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Fingerprints for spin-selection rules in the interaction dynamics of O2 at Al(lll) 

Christian Carbogno,^ Jorg Behler,^ Axel GroB/ and Karsten Renter^ 

' Institut fiir Theoretische Chemie, Universitdt Ulm, D-89069 Ulm, Germany 

Lehrstuhl fiir Theoretische Chemie, Ruhr-Universitat Bochum, D-44780 Bochum, Germany 

^ Fritz- Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany 

(Received 18th June 2008) 

We performed mixed quantum-classical molecular dynamics simulations based on first-principles 
potential-energy surfaces to demonstrate that the scattering of a beam of singlet O2 molecules at 
Al(lll) will enable an unambiguous assessment of the role of spin-selection rules for the adsorption 
dynamics. At thermal energies we predict a sticking probability that is substantially less than unity, 
with the repelled molecules exhibiting characteristic kinetic, vibrational and rotational signatures 
arising from the non-adiabatic spin transition. 
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Chemical interactions are governed by various selec- 
tion rules, among which the ones due to Wigner formu- 
late the constraint of overall spin-conservation 1]. For 
the ubiquitous oxygen gas-phase chemistry a well-known 
implication of these spin-selection rules is the pronounced 
inertness of the O2 molecule in its triplet ground state, if 
the other reactant and the product are spin singlets. In- 
triguingly, the dissociative adsorption of oxygen at metal 
surfaces represents exactly such a spin constellation, but 
there spin-selection rules are rarely perceived as poten- 
tial reason for a low reactivity. A notable exception is 
the long-time enigmatic low initial sticking probability 
of thermal O2 molecules at Al(lll) [3, |3|, which was 
recently explained in terms of a corresponding strongly 
non-adiabatic dissociation dynamics 4] : Due to the hin- 
dered spin transition impinging O2 molecules maintain 
their initial triplet configuration and are already repelled 
back into the gas-phase at rather large distances from the 
surface where other mechanisms such as charge transfer 
from the metal are still quite weak. This non-adiabatic 
picture of the dissociation process was substantiated with 
extensive locally constrained density-functional theory 
(DFT) [5|, la] calculations, which allow to localize the 
triplet spin at the oxygen molecule and thereby give ac- 
cess to the corresponding spin-triplet potential-energy 
surface (PES). Ensuing molecular dynamics (MD) simu- 
lations confining the trajectories of the approaching O2 
molecules to this spin-triplet PES revealed indeed a re- 
duced sticking probability So{E) at low kinetic energies 
E, in semi-quantitative agreement with the experimental 
data and in gross contradiction to the constant and essen- 
tially unit sticking probability obtained from equivalent 
MD simulations on the adiabatic PES.yi 

Despite this significant progress to firmly establish the 
O2 scattering at Al(lll) as a first benchmark example for 
the importance of spin selection rules at metal surfaces, 
the limitations due to the semi-local DFT exchange- 
correlation functionals available for the mapping of the 
corresponding high-dimensional PESs require further at- 
tention. At present, we can not rule out that an improved 



description of electronic exchange and correlation would 
e.g. not yield barriers on the adiabatic PES and there- 
with a similarly reduced sticking probability within an 
entirely adiabatic framework. While this obviously dic- 
tates further research along these lines \m, LZ| , one also has 
to recognize that the hitherto targeted triplet O2 stick- 
ing probability is not necessarily an observable that is 
most sensitive to the non-adiabatic spin flips themselves. 
In fact, due to the position of the triplet PES barriers 
at rather large distances from the surface, the shape of 
So{E) is to a large extent governed by the motion on 
this PES alone, and it is for this reason that the previous 
modeling could do even without explicitly accounting for 
the actual spin transitions themselves. (J] 

In this Letter we therefore suggest to change the focus 
to a complementary observable that depends sensitively 
on the spin coupling and which will constitute unequivo- 
cal and measurable fingerprints for the existence or non- 
existence of spin-selection in this system: the initial stick- 
ing probability for spin-singlet O2 molecules, S*{E). As 
we will demonstrate on the basis of first-principles calcu- 
lations, this probability is significantly lower than unity 
at low kinetic energies due to the hindered, but non-zero 
transitions between the different spin configurations of 
the impinging O2 molecules. In contrast to the case of the 
spin-triplet sticking, the repelled molecules leave the sur- 
face this time in the flipped spin-state, i.e. they impinge 
as spin singlets, but leave as spin triplets. The important 
twist is thus that the repelled molecules have gained the 
triplet-singlet energy, which in the experiment is about 
AEts ~ 1 eV. From the calculated redistribution of this 
sizable excess energy into translational, rotational and 
vibrational degrees of freedom we predict that these ex- 
iting molecules will therefore lead to signatures that will 
clearly allow to discriminate them as fingerprints for the 
non-adiabatic dynamics, even if the molecular beam em- 
ployed does not entirely consist of singlet O2 molecules. 

Our first-principles calculations of the sticking coeffi- 
cients are based on the general "divide and conquer" - 
style approach and on exactly the same computational 
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FIG. 1: Calculated singlet- and triplet potential energy along 
the minimum energy path for a triplet O2 molecule impinging 
side-on over the Al(lll) fee hollow site, cf. the schematic top 
view shown as inset. The energy zero corresponds to a free 
O2 molecule in its triplet ground state. At corresponding dis- 
tances Z — » cx) far away from the Al(lll) surface, the singlet 
energy is higher by A_Ets ~ 1 eV. Additionally shown is the 
spin coupling matrix element Vc along the path. Obviously, 
the transitions due to this coupling are strongest at the cross- 
ing seam of the singlet and triplet PESs, which is located at 
rather large distances Z > 2 A away from the surface. 



setup for the 02/Al(lll) system as detailed in Ref. [^. 
In short, this comprises an extensive mapping of the sin- 
glet and triplet PESs along the six molecular degrees of 
freedom using DFT within the generalized gradient ap- 
proximation (GGA) [S]], followed by an interpolation us- 
ing neural networks [9|, llO| in order to obtain the con- 
tinuous representation required for the MD simulations. 
Locally constrained DFT calculations enforcing a triplet 
spin at the oxygen molecule are employed to map the 
triplet PES (corresponding to the ^E~ ground state at 
large distances from the surface) 6]. Due to the improper 
multiplet representation given by DFT with jellium based 
exchange-correlation functionals, the same approach but 
enforcing a singlet spin at the oxygen molecule can not be 
used to obtain an appropriate description of the lowest- 
lying singlet PES (corresponding to the ^ Ag excited state 
at large distances from the surface). As discussed be- 
fore [Jl an energetically reasonable approximation to this 
singlet PES is instead obtained through spin-unpolarized 
DFT calculations, yielding a triplet-singlet energy for the 
free O2 molecule of AE'ts — 1-2 eV. Despite this approxi- 
mate nature of the singlet state, the different spinors still 
ensure the orthogonality of the two electronic configura- 
tions, which can thus be viewed as spanning a kind of 
two-dimensional diabatic basis. 

Figure [T] illustrates the two PESs relevant for the O2 
sticking along a one-dimensional cut through configura- 
tion space, corresponding to the minimum energy path 
for a triplet O2 molecule impinging side-on over the 
Al(lll) fee hollow site. Whereas the preceding work con- 



sidered only motion confined to one of these electronic 
states, we now explicitly account for an excited state 
dynamics and non-adiabatic transitions by performing 
mixed quantum-classical MD simulations based on the 
fewest switches algorithm pjj, [l^ . In this surface hop- 
ping framework, the nuclear degrees of freedom are inte- 
grated classically on one PES at each time step. Simul- 
taneously the density matrix including both electronic 
states is calculated by integrating the time-dependent 
Schrodinger equation along this trajectory. Transitions 
from one PES to another are introduced in such a way 
that for a large number of trajectories the occupation 
probabilities given by the density matrix are achieved 
within the smallest number of switches possible. The ac- 
tual spin coupling matrix element Vc at each point of the 
configuration space spanned by the six molecular degrees 
of freedom is calculated by taking advantage of the also 
available adiabatic spin-polarized ground-state PES [J]. 
Inverting the diagonalization of the diabatic Hamiltonian 



H 



(Et Vc 

\Vc* Es 



where Et and Eg are the triplet and singlet energies re- 
spectively, then uniquely determines the absolute value of 
the matrix elements Vc in all six considered dimensions, 
cf. Fig. [TJ For the present two-state problem, this thus 
enables a computationally efficient determination of the 
non-adiabatic couplings without explicitly requiring to 
evaluate the underlying wave function dynamics [l3| . We 
note, however, that particularly because of the approx- 
imate nature of the singlet PES, this procedure likely 
overestimates the Vc and we will critically discuss this 
point below. 

As a natural first target we use this approach to recom- 
pute the sticking curve for impinging triplet O2 molecules 
that was the central quantity addressed in the preced- 
ing study [4|. More than 2000 trajectories with random 
initial molecular orientations were calculated for each ki- 
netic energy, in order to get a reliable statistics for this 
averaged quantity. As suspected, we obtain a result that 
is very similar to the curve published in Ref. [J], i.e. 
a So{E) that is close to zero for thermal molecules, in- 
creases monotonically with increasing kinetic energy and 
approaches 100% for E > 0.5 eV. While this reconfirms 
that hindered spin transitions lead indeed to a lowering 
of So{E) at low kinetic energies as observed experimen- 
tally, we also recognize that the explicit consideration of 
spin transitions in our mixed quantum-classical dynamics 
barely modifies the result obtained previously with a dy- 
namics that was entirely confined to the triplet PES [J]. 
As apparent from Fig. [U the lowest-energy molecules are 
repelled by the barriers on the triplet PES already at dis- 
tances from the surface, where transitions to the singlet 
state are energetically forbidden. This determines the 
characteristic "S" -shape of the sticking curve, and the fi- 
nite spin transitions only lead to a small increase of So {E) 
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FIG. 2: Initial sticking coefficient at normal incidence for 
singlet oxygen molecules, S^iE), as computed for a dynam- 
ics confined to the singlet PES (V) and by accounting for 
non-adiabatic transitions (D). The dashed line represents the 
sticking coefficient for triplet oxygen, simply shifted by A_Ets 
(see text). 



at intermediate energies between those of the triplet- 
singlet crossing seam and the triplet barriers, marked as 
-S'cross a-nd -Ebarr i^ Fig- HI The initial sticking coefficient 
for triplet O2 molecules at low kinetic energies is thus 
not very sensitive to the details of the non-adiabatic cou- 
pling. This is aggravated by the fact that in this energy 
window the fraction of repelled molecules, (1 — So{E)), 
consists entirely of O2 molecules in their triplet ground- 
state. Imagine that the absence of barriers on the hith- 
erto calculated adiabatic PES was an artifact of the em- 
ployed semi-local DFT exchange-correlation functional, 
and the real adiabatic PES exhibited barriers just of the 
same order as the present triplet PES. This would equally 
yield an "S" -shaped So{E) curve, explaining the experi- 
mental data in an entirely adiabatic framework. Further- 
more, since also here the repelled molecules are in their 
spin-triplet state, it would be impossible to distinguish 
whether they had been repelled adiabatically or as the 
result of a hindered spin-transition. 

This situation is markedly different, if we move on to 
the sticking curve for a beam of singlet O2 molecules. 
From the purely attractive character of the corresponding 
singlet PES, cf. Fig. [H one would intuitively expect a 
sticking curve of S*{E) « 1 at all kinetic energies, and 
this is indeed what is obtained when confining the MD 
simulations to a motion on the singlet PES as shown in 
Fig. [51 Surprisingly, accounting for the finite singlet- 
triplet coupling leads now to a significant reduction of 
the sticking probability at low kinetic energies, cf. Fig.[2l 
This is due to the fact that some of the impinging singlet 
molecules suffer a spin transition near the crossing seam 
of the singlet and triplet PES, and are back-scattered 
on the triplet PES, which indeed exhibits some barriers 
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FIG. 3: Percentage of the electronic potential energy A_Bts 
redistributed into translational (o), rotational (D) and vibra- 
tional (A) energy of the refiected O2 molecules. Note that all 
these molecules leave the surface as spin-triplets. 



larger than 1 eV. At increasing kinetic energies, this effect 
is reduced as more molecules can overcome these effective 
barriers, but also because the singlet-triplet transition 
probability is reduced at higher velocities 14 1 . 

For the lowest kinetic energies, however, this transi- 
tion probability is almost 100%, i.e. virtually all im- 
pinging singlet molecules cross to triplet when approach- 
ing the singlet-triplet PES crossing seam. One might 
therefore expect that at these energies the singlet stick- 
ing curve simply resembles the triplet sticking curve, 
just shifted by the vacuum triplet-singlet energy, i.e. 
S*{E) w So{E + AEts)- Again surprisingly, this is not 
the case and the surface hopping simulations yield an 
even more pronounced reduction of S* (E) at the lowest 
kinetic energies, cf. Fig. [21 This is rationalized by re- 
alizing that the coupling matrix elements are not only 
a function of the molecule-surface separation, but also 
of the other molecular degrees of freedom. As a conse- 
quence, the excess energy released upon transition from 
singlet to triplet PES is redistributed onto these differ- 
ent degrees of freedom according to the direction of the 
non-adiabatic coupling vector [llj . Correspondingly, not 
the entire potential energy stored in the electronic sin- 
glet configuration is available for the propagation along 
the reaction path and the dissociation probability is even 
further reduced. 

A fingerprint for this energy redistribution during the 
non-adiabatic singlet-triplet transition is given by the 
amount of energy stored into the translational, vibra- 
tional and rotational degrees of freedom of the reflected 
molecules. Figure [3l summarizes this data, and for this 
it is important to realize that we find essentially all re- 
flected molecules to leave the surface in the triplet spin 
state. Due to the redistributed excess energy AEts, the 
reflected molecules become translationally, vibrationally 
and rotationally hot, i.e., upon reflection they gain a sig- 



nificant amount of energy in all molecular degrees of free- 
dom. Their energy distribution thus provides unequivo- 
cal signatures for the non-adiabatic spin transition. As 
apparent from Fig. [31 this holds in particular for ther- 
mal molecules, for which roughly one third of Ai^xSj 
thus more than 0.3 eV, is stored in the translational mo- 
tion [iSJ I . Even if the experimental molecular beam does 
not completely consist of singlets, it would therefore eas- 
ily be possible to discriminate non-adiabatically reflected 
molecules by their kinetic energies alone. 

In contrast to the hitherto much studied sticking of 
triplet O2 molecules we thus find the scattering of sin- 
glet molecules at Al(lll) to yield observables in form of 
the reflected molecules that would clearly flag a poten- 
tially non-adiabatic dynamics governed by spin selection 
rules. If spin transitions were not hindered at all close 
to the surface, essentially all impinging singlet molecules 
would stick, since they would quickly release the potential 
energy stored in the electronic configuration and there- 
with be able to overcome even modest barriers on the 
adiabatic PES. With selection rules operating, our cal- 
culations predict instead that for example ~ 20% of an 
impinging beam of thermal singlet molecules will be re- 
pelled with characteristic translational, vibrational and 
rotational signatures. In this respect, we come back to 
the point that our coupling matrix elements likely over- 
estimate the real transition probabilities, most notably 
because of the approximate singlet PES and the dia- 
batic Hamiltonian inversion procedure involving the fully 
adiabatic PES. The resulting too strong destabihzation 
of the singlet-state is particularly pronounced at large 
molecular distances from the surface due to the spurious 
charge transfer into the O2 molecule contained in the sin- 
glet PES [4|. However, even if we reduce the Vc to the 
value of the vacuum spin-orbit coupling that then yields 
a gas-phase singlet lifetime that is orders of magnitude 
longer than the duration of the entire surface scatter- 
ing event, and consider the absolutely unrealistic case 
that the Vc stays constant at this minute value even at 
any closer distance to the surface, we still obtain a re- 
duced, but nevertheless clearly measurable percentage of 
reflected molecules with distinct kinetic, vibrational and 
rotational properties [16[. Moreover, the redistribution 
of the excess energy into the different degrees of freedom 
and in particular into the translational motion of the re- 
pelled molecules is virtually not changed, so that our 
conclusions in terms of unequivocal fingerprints remain 
completely untouched. 

Summarizing, we have performed mixed quantum- 
classical dynamical simulations based on first-principles 
potential-energy surfaces to demonstrate that the scat- 
tering of a beam of singlet O2 molecules at Al(lll) will 
enable an unambiguous assessment of the role of spin- 
selection rules for the adsorption dynamics. We predict 
a sticking probability at low kinetic energies that is sub- 
stantially less than unity, with the repelled molecules 



exhibiting characteristic kinetic, vibrational and rota- 
tional signatures of a non-adiabatic spin transition that 
will allow to discriminate them even if the molecular 
beam does not purely consist of spin-singlets. Such ex- 
periments are well feasible as recently demonstrated by 
Burgert et al. [17|], who already pointed out the impor- 
tance of spin-selection rules for the interaction of singlet 
oxygen interacting with small anionic clusters formed of 
~ 10 — 20 Al atoms. With the here predicted fingerprints, 
it thus remains to future experimental observation to de- 
cide whether spin-selection rules are also operative at the 
extended metallic surface or not. If yes, the present ap- 
proach (mixed quantum-classical simulations based on 
PESs from constrained DFT) sets the conceptual stage 
to address corresponding effects in the reaction dynamics 
of more complex systems. 
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